Relationship between low and high frequencies in Delta Scuti 
stars: Photometric Kepler and spectroscopic analyses of the rapid 

rotator KIC 8054146. 1 
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ABSTRACT 

Two years of Kepler data of KIC 8054146 ((5Sct/7Dor hybrid) revealed 349 statistically sig- 
nificant frequencies between 0.54 and 191.36 cd -1 (6.3 /iHz to 2.21 mHz). The 117 low frequencies 
cluster in specific frequency bands, but do not show the equidistant period spacings predicted for 
gravity modes of successive radial order, n, and reported for at least one other hybrid pulsator. 
The four dominant low frequencies in the 2.8 to 3.0cd _1 (32 to 35 /iHz) range show strong ampli- 
tude variability with timescales of months and years. These four low frequencies also determine 
the spacing of the higher frequencies in and beyond the 6 Set pressure-mode frequency domain. 
In fact, most of the higher frequencies belong to one of three families with spacings linked to 
a specific dominant low frequency. In the Fourier spectrum, these family regularities show up 
as triplets, high-frequency sequences with absolutely equidistant frequency spacings, side lobes 
(amplitude modulations) and other regularities in frequency spacings. Furthermore, within two 
families the amplitude variations between the low and high frequencies are related. We conclude 
that the low frequencies (gravity modes, rotation) and observed high frequencies (mostly pres- 
sure modes) are physically connected. This unusual behavior may be related to the very rapid 
rotation of the star: from a combination of high and low-resolution spectroscopy we determined 
that KIC 8054146 is a very fast rotator (v sin i = 300 ± 20 km s" 1 ) with an effective temperature 
of 7600 ± 200 K and a surface gravity log g of 3.9 ± 0.3. Several astrophysical ideas explaining the 
origin of the relationship between the low and high frequencies are explored. 

Subject headings: Stars: oscillations (including pulsations) - Stars: rotation - Stars: individual: 
KIC 8054146 - Stars: variables: delta Scuti - Stars: abundances 

1. Introduction 
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Kepler observations have completely changed 
our understanding of stellar pulsations in the S Set 
instability strip. Earlier ground-based data had 
shown that 7 Dor stars, which pulsate in frequen- 
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cies typically below 5d _1 , lie in a fairly small re- 
gion on, or just above, the main sequence that 
partly overlaps the cool edge of the 8 Set insta- 
bility strip. A few stars were known in which 
both 7 Dor and 8 Set pulsations were present (the 
so-called hybrids). Kepler photometry has shown 
that practically all A and F stars, whether they are 



5 Set stars or not, have low frequencies (jBalona 
20111 ). Thus the term 'hybrid' has lost its mean- 
ing since nearly all 8 Set stars are hybrids. The 
convective blocking mechanism, which is thought 
to drive 7 Dor pulsations, ceases to be effective for 
stars hotter than about 7500 K, yet low frequen- 
cies are found in the hottest A-type stars. It is, 
therefore, not possible to fully understand the low 
frequencies in the A-type stars as pulsation. 

Statistical analyses of the low frequencies 
in A-type stars show that the frequencies are 
compatibl e with the expected rotational fre- 
quencies (jBalona 1 120111) . In fact, for the 8 Set 
star KIC 9700322, a secure identification of low- 
frequency peaks with stellar rotation was estab- 
lished from a variety of techniques: stellar line 
width from spectroscopy, the spacing of an ob- 
served £=2 p-modc quintuplet and the time scales 
of the small amplitude modulation of the dom - 



inant modes (jBreger et al. l201lt i reger 1 1201 If 



The low and high frequencies in hybrid pulsators 
are also related through additional physical mech- 
anisms, which are presently not fully understood. 
This relationship is seen in a number of 6 Set stars 
measured by space telescopes. An excellent exam- 
ple is the C0R0T star ID 105733033 studied by 
Chapellier et al.l ([20121 ). which shows a number of 
couplings between the low and high frequencies. 
Another star showing a number of numerical re- 
lationships between low and high frequencies is 
KIC 8054146, the star examined in the present 
paper. This star is interesting from another point 
of view as well: its very high rotational velocity 
(see below). 

In general, astcroseismic studies of rapidly ro- 
tating stars of spectral type A and early F have 
been less successful than those of slow rotators. 
The reason is both theoretical and observational: 
rapid rotation cannot yet be reliably treated in 
theoretical modeling, while spectroscopically the 
extreme width of the spectral lines of the rapidly 
rotating stars makes detailed abundance anal- 
yses very difficult. Furthermore, observational 



studies (e.g., iBreger l2000h have shown that in 
8 Set stars low rotation is a prerequisite for high- 
amplitude radial pulsation, while the rapidly ro- 
tating 8 Set stars have small photometric ampli- 
tudes with mainly nonradial pulsation. 

In this paper we shall study the asteroseismic 
properties of a very rapidly rotating A/F star, 
KIC 8054146, which is also a hy brid pulsator with 
both 8 Set and 7 Dor properties (jUvtterhoeven et al 
20111 ). i.e., a star with reported gravity and pres- 
sure modes. 

2. Kepler observations of KIC 8054146 

The Kepler mission is designed to d etect Earth- 
like p lanets around solar- type stars ( Koch et al.l 
20101 ). To achieve that goal, Kepler is continuously 



monitoring the brightness of over 150000 stars for 
at least 6yr in a 115 square degree fixed field of 
view. Photometric results show that after one 
year of almost continuous observations, pulsation 
amplitudes of 2 ppm (parts-per-million) are eas- 
ily detected in the periodogram for stars brighter 
than V = 11 mag. Two modes of observation are 
available: long cadence (LC, 29.4-min exposures) 
and short-cadence (SC, 1-min exposures). With 
short-cadence exposures it is possible to observe 
the whole frequency range from to 200 cd -1 seen 
in 8 Set, 7 Dor, and hybrid stars pulsating with 
both the 8 Set and 7 Dor frequencies. 

KIC 8054146 (K p = 11.33) was extensively ob- 
served with the Kepler satellite in short cadence 
during quarters Q2.3 (BJD 245 5064 - 245 5091) as 
well as Q5 to Q10 (BJD 245 5276 - 245 5833), cov- 
ering a time span of 27 d and 557 d, respectively. 
The large amount of accurate data covering more 
than a year leads to an unprecedented frequency 
resolution, permitting us to precisely examine the 
different frequency patterns (separations). The 
coverage is ideal to examine the relation between 
the low-frequency variability (possibly caused by 
gravity modes and stellar rotation) and the higher- 
frequency (5 Set pressure modes. 

The Kepler short-cadence photometry with a 
measurement per minute is available as calibrated 
and uncalibrated data. To avoid the known sys- 
tematic errors present in the calibrated data, we 
only used the uncalibrated data. We performed 
our own reductions eliminating the zero-point 
trends, outliers, jumps and known spurious ef- 
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Fig. 1. — Typical light curve covering 0.75 days out of a total of more than 580 days. The plot demonstrates 
the presence of both low- and high-frequency groups. The average deviation per single point from the fit is 
0.18 ppt. 



fects such as those caused by the reaction wheels 
(spurious frequencies near 0.3 cd -1 ). Our inde- 
pendently develope d reduction p r ocess is similar 
to that published bv lGarcia et al. (l201ll ) for solar- 
type stars and the description will not be repeated 
here. Despite the careful reductions, any detec- 
tions of very low frequencies below 0.3 cd _1 should 
be treated with caution, but this does not affect 
the present conclusions. For Q8, a known spu- 
rious frequency, drifting from 31.5 to 32.2 cd -1 
with amplitudes around 60 ppm was taken out on 
a daily basis. 

The short-cadence data were supplemented by 
long-cadence data from quarters Ql to Q4 (BJD 
245 4964 - 245 5275). These measurements of two 
points per hour were reduced by us in a similar 
manner to the short-cadence data. A compari- 
son of the resulting LC and SC frequency spec- 
tra for the important Q2.3 time period gave an 
excellent average agreement to ± 3 ppm for the 
low-frequency peaks, corresponding to the size of 
the formal uncertainties. Visual comparisons for 
subsequent quarters with both types of data con- 
firmed that long-cadence data can indeed be used 
for studies of low-frequency amplitude variability 
(only). This result is in agreement with a discus- 
sion of the cha racteristics of the Kepler data by 
Murphy I ([20121 ). Nevertheless, we emphasize that 
in this paper the LC data are used only to exam- 
ine the low-frequency amplitude variability from 



Ql to Q4. 

3. Frequency analyses 

The Kepler short-cadence data of KIC 8054146 
were analyzed with th e statistical package PERI0D04 
(jLenz fc Bregerl 120051 ) . This package carries out 
multifrcquency analyses with Fourier as well as 
least-squares algorithms and the determination of 
significance of a peak does not rely on the assump- 
tion of white noise. 

Following the standard procedures for ex- 
amining the peaks with PERI0D04, we have 
ado pted the amplitud e signal/noise criterion of 
4.0 (|Breeer et al.lll993l) . This standard technique 



is modified for all our analyses of accurate satel- 
lite photometry: the noise is calculated from 
prewhitened data because of the huge range in 
amplitudes of three orders of magnitudes. The 
extensive data led to very low root-mean-squarc 
noise levels in the Fourier spectrum: lppm (0.5 - 
7cd~ 1 range), 0.75 ppm (7-20cd _1 ), 0.6ppm (20 
- SOcd" 1 ), 0.55 ppm (30 - eOcd" 1 ), and 0.5ppm 
(60 - 200 cd- 1 ). 

We undertook extensive computations to elimi- 
nate the effects of amplitude and (small) frequency 
variations, while preserving the extremely high 
frequency resolution. A large number of tests in- 
volving various amounts of prewhitening of the 
dominant frequencies were also undertaken with 
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Fig. 2. — The amplitude spectra of KIC 8054146. The top panel shows the spectrum of the data covering 
two years. To eliminate the effects of the frequency drifts in the star, the bottom panel uses Fourier 
spectra covering three months only. Furthermore, in order to demonstrate the many frequencies with small 
amplitudes, the panel shows amplitudes logarithmically. Note the large frequency range over which the star 
is variable. 



thousands of hours of computations. 

We detected 349 statistically significant fre- 
quencies between 0.54 and 191.36 cd -1 . There is a 
clear separation in the density of frequency peaks 
between low and high frequencies: 116 frequencies 
have values smaller than 6cd _1 , there is a single 
peak at 6.71 cd -1 , while 232 peaks are found be- 
tween 7 and 192 cd -1 (the high-frequency group). 
The results are shown in Figure [2] which presents 
the Fourier spectra in two different ways. The 
top panel ignores the small period and large am- 
plitude variability of all frequencies and presents 
the average solution. The lower panel represents 
the average of separate Fourier analyses for each 
quarter and depicts the result with logarithmic 
amplitudes in order to emphasize the many low- 



amplitude peaks. 

The frequencies with average amplitudes of 
5.0 ppm or larger are listed in Tabled] Due to the 
strong amplitude variability, these average ampli- 
tudes were not computed from a single solution, 
but were obtained from averaging the amplitudes 
of the three Q2, Q5-7, and Q8-10 time periods. 
The Q2 data cover only 27 days, but were given 
equal weight for determining the averages because 
of the dominance of fi and all its related sequences 
in the early quarters only. On the other hand, the 
lower frequency resolution of the Q2 data forced us 
to omit some close frequency pairs from the Q2 so- 
lution, explaining a small number of blank spaces 
in the table. The 120 frequencies in the table can 
only be a subset of the statistically significant frc- 
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qucncies and by themselves cannot reveal all the 
rich regularities present in the star. These will be 
covered below. 

Table Q] also lists a number of identifications 
such as combination frequencies, harmonics and 
other sequences, which will be covered below. We 
have restricted the listed identifications to only 
the main patterns, mostly involving the dominant 
frequencies fi to fi6- Multiple listings for each 
frequency component in a combination were also 
avoided in the interest of clarity. We also note 
the excellent frequency accuracy: the predicted 
and observed frequencies of the combinations fit 
within a few 0.0001 cd" 1 . 

Are there hundreds or thousands of additional 
frequencies which we have not detected? If we ex- 
clude the low-frequency region, in which the noise 
is definitely not white, we consider our amplitude 
limit of 2 to 3 ppm quite restrictive and can rule 
out a large number of additional modes with pho- 
tometric amplitudes of this size or larger. In addi- 
tion, no statistically significant peaks at higher fre- 
quencies up to the Nyquist frequency at 368 cd -1 
were detected. 
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Fig. 3. — The distribution of the low-frequency 
peaks in period (= 1/frequency) . Period, rather 
than frequency, is used here since gravity modes 
are expected to be equally spaced in period, rather 
than frequency. Equidistant period spacings have 
not been detected. Note that the peaks cluster 
mainly in four specific regions. 



4. The low frequencies 

The 117 low frequencies detected by us arc not 
randomly distributed. They are mostly clustered 



in specific frequency bands: 1.6 - 1.8, 2.6 - 3.0, 
3.5 - 3.7, and 5.2 - 6.0 cd -1 . We note that the 
last two bands are located at approximately twice 
the frequencies of the first two bands, although 
only some individual peaks obey a strict f, 2f rela- 
tion. The clustering of the peaks is demonstrated 
in Figure[3l Since gravity modes of adjacent radial 
order, n, are expected to show equal spacing in pe- 
riod in the asymptotic limit, Figure [3] shows the 
distribution in period (inverse frequency). A cur- 
sory examination already reveals that equidistant 
period (or frequency) spacing is not dominant. 

For the hybrid pulsa tor, CoRoT ID 105733033, 
Chapellier et al.l ( 2012t ) have discovered a remark- 
able sequence of asymptotically spaced frequen- 
cies with an equidistant period spacing of 0.0307 
days. They interpret this to be caused by £ = 2 
gravity modes of different radial order, n. Since 
for KIC 8054146 we also see a large number of 
closely spaced peaks at low frequencies, the ques- 
tion arises whether such gravity modes are also 
present. In KIC 8054146 the amplitudes of the 
low-frequency peaks are a factor of 10 - 100 
smaller, but the extensive data set makes such a 
search possible. 

However, a more cautious approach might be in 
order. Some aliasing exists in our Kepler data due 
to small gaps in coverage. Furthermore, the slow 
amplitude variability in this star can also lead to 
close frequencies. These effects need to be mini- 
mized for the the low-frequency region with many 
close frequency peaks. 

Consequently, we have repeated our previous 
analysis with some refinements: To reduce any 
effects of aliasing caused by the large gap be- 
tween Q2 and Q5 as well as the small gaps in 
the Q5 to Q10 data, we omitted the 27d of the 
Q2 data. This reduced the aliasing to four small 
peaks separated from the main frequencies by less 
than 0.01 cd _1 . To reduce the effects of amplitude 
variability of the dominant peaks, we calculated 
multifrcquency solutions for the 10 dominant low 
frequencies with PERI0D04 allowing for quarterly 
amplitude and phase changes. The prewhitened 
data were then used to find additional, statisti- 
cally significant peaks with amplitudes of 5 ppm 
or larger. 

We then examined the four regions for regular 
frequency as well as period spacings of the individ- 
ual peaks by constructing histograms of the ob- 
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served frequency and period differences. Again, 
regular patterns could not be detected, in par- 
ticula r, the spacings observed by IChapellier et al 



(|2012h in another hybrid star were not found. 
We conclude that for KIC 8054146, gravity modes 
with equidistant period spacings were not de- 
tected. 

The four dominant low frequencies in the 2.8 
to 3.0cd _1 (32 to 35/iHz) show strong amplitude 
variability with time scales of months and years 
(Figurc[4}. An example is the peak at 2.814cd~ 1 
(fi), which reaches its maximum value during Q3, 
slowly decreases its amplitude during Q5 and is 
almost absent from Q6 to Q9. The 2f term fol- 
lows this behavior. The 2.98 cd _1 peak, however, 
appears in Q3, slowly grows and starts to disap- 
pear during Q9 and Q10. The behavior of the 
2.93/5.86cd~ 1 group of peaks is more complex, 
and will be discussed below together with its re- 
lated family of high-frequency modes. 
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Table 1: Multifrcqucncy solution of KIC 8054146, patterns and identifications 
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T: triplets, S: equidistant picket-fence series. The table contains only frequencies with amplitudes of > 5.0 ppm. 
Almost all amplitudes are variable and the averages are listed. Formal errors are 2 ppm or smaller. 
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Fig. 4. — Amplitude variability of the dominant 
low frequencies over 2.5 years. The dominant 
peaks in 5 - 6cd _1 range are mostly multiples 
and combinations of the 2.5 - 3.0 cd -1 peaks. 
The amplitude variations of the 2f multiples are 
not in step with the amplitude variations of f. 
Note also the strong amplitude decrease of the fi 
(2.814 cd -1 ) peak after Q5 and the "return" at 
2.806 cd" 1 after Q8. 



5. The high frequencies 

We have detected 232 statistically significant 
frequency peaks between 7 and 191 cd -1 . This 
large range in frequency is unusual for S Set stars. 
Not surprisingly, some of the frequency peaks are 
expected linear combinations of high frequencies 
and harmonic terms, but cannot explain most of 
the high frequencies. 

The dominant high frequencies are three groups 
of triplets (Ti to T 3 ), the difference between the 



Ti and T 2 triplets (= 25.952cd" 1 ), and a 'lonely' 
peak at 17.349 cd -1 . The triplets show very strong 
amplitude variability with a timescale of a year or 
longer. The most striking property of the detected 
frequency peaks is their separation in frequency, 
which often is identical to the previously detected 
low-frequency peaks. 



Frequency spacing (^iHz) 
20 30 40 50 60 




3 4 5 

Frequency spacing (c/d) 



Fig. 5. — Histogram of the frequency spacings 
among the high frequencies (7 to 200 cd -1 ) using 
a narrow bin size of 0.0033 cd -1 . The most com- 
mon spacings are found to be numerically equal 
to the dominant low frequencies, which have been 
marked in the diagram. The diagram shows that 
the low frequencies in the g-modc region strongly 
affect the spacings of a large number of high fre- 
quencies in the p-modc domain. 

A powerful method to look for regularities in 
frequency spacing consists of forming all possible 
frequency differences between all significant fre- 
quency peaks and to examine the histograms. The 
huge number of computed differences has to lead 
to some accidental agreements. In the present 
case, the high accuracy of the detected frequencies 
permits us to choose extremely narrow histogram 
boxes. Furthermore, we restrict the analysis to 
spacings smaller than 7cd _1 . We have calculated 
all the frequency spacings between the 232 high- 
frequency peaks between 7 and 200 cd -1 . A his- 
togram was computed with a very small box size 
of 0.0033 cd^ 1 . The resulting histogram with 2100 
boxes is shown in Figure [5] We regard the result 
as remarkable. 

The four dominant low frequencies, fi to fj, 
account for 136 frequency spacings among the 
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232 high-frequency peaks. A few of these spac- 
ings (within the narrow adopted range of only 
0.0033cd _1 ) may be accidental. However, we be- 
lieve the real number of exact spacings correspond- 
ing to observed low frequencies to be even higher 
because of our restriction to only simple agree- 
ments, e.g., spacings involving the (fi+f/i) sum are 
ignored by our technique. 

The preferred spacings involving f'i to £4 can be 
easily understood in terms of three families of pul- 
sation modes. Each family behaves differently and 
may require different astrophysical explanations. 

5.1. The f2, f"3 family and triplets 

Figure |6] shows the two triplets, Ti and T2, and 
their relationship to three low frequencies and the 
25.9517cd _1 peak. All shown peaks have high am- 
plitudes with a corresponding high accuracy of the 
frequency values (± 0.0002 cd _1 ). While data cov- 
ering only one quarter would suggest equal sep- 
arations, the two-year data made it possible to 
detect f2 and f3, which are clearly separated by 
0.0040 ± 0.0003 cd _1 . This small difference has 
strong implications. We note: 

Frequency in |.cHz 
100 200 300 400 500 600 700 800 
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2.9299 
<f 2 ) 



10 



20 



30 40 50 
Frequency in c/d 



60 



70 



Fig. 6. — The spacings of the two high-frequency 
triplets and their relation to dominant low fre- 
quencies. The two triplets are separated by 
25.951.7cd~ 1 , which is also detected as a sepa- 
rate frequency peak. Also, the spacings between 
the triplets correspond to a close low-frequency 
pair and the sum of that pair. The separate Ti 
and T2 triplet components share similar amplitude 
changes. 



(i) The two triplets have nearly identical prop- 
erties, including amplitude variability. Their fre- 
quency separation of 25.9517cd _1 is also seen as 
a strong separate peak. 

(ii) The frequency differences between the cen- 
tral and right peaks correspond to an observed, 
dominant low frequency peak at 2.9300 cd _1 (fij). 

(iii) The frequency differences between the cen- 
tral and left peaks correspond to an observed, 
dominant low frequency peak at 2.9340 cd _1 (fa). 

(iv) The frequency differences between the left 
and right peaks correspond to an observed, dom- 
inant low frequency peak at 5.8641 cd _1 (f2+fa). 
This peak is not a harmonic of cither f2 or f3. 

(v) The amplitude changes of the f2, f3, and 
(f2+fs) peaks mirror the amplitude changes of 
the individual components of the triplets. The 
decrease of the f2 and fj strengths after Q7 are 
closely mirrored by the decrease of the amplitudes 
of the central triplets components. On the other 
hand, the increase of the (fi+f2) amplitude after 
Q7 is mirrored by the increase in the non-central 
components of the triplets. (The lengthy detailed 
analysis of this phenomenon, which includes anal- 
yses of weekly amplitude and phase changes, will 
be the subject of a later report.) 

The observed measurements are incompatible 
with an interpretation that the 5.86 cd -1 peak 
is a 2f harmonic of a mode near 2.93 cd _1 and 
may challenge such interpretations in other hybrid 
stars. In Section [7] we will explore the possibility 
that the triplets are the counterparts of the I = 1 
triplets commonly seen in slowly rotating stars and 
the possibility of combination modes. 

5.2. Equidistant S sequences involving fi 

Between 76 and 102cd _1 , the bottom panel 
of Figure [5] shows five equidistantly spaced peaks 
with additional smaller peaks exactly halfway be- 
tween them. The separation is exactly equal to the 
dominant low frequency, fi, at 2.814cd~ 1 . The 
frequencies of the sequence are not exact multi- 
ples of fi , but obey a relation of f = frj + n fi , where 
f = 2.852, fi = 2.8142, and n ranges from 25 to 34. 
(Note that the high precision of the spacing, fx, is 
achieved by dividing the detected high frequencies 
by n. Furthermore, fo is the mathematical inter- 
cept of the relation and not a detected peak.) We 
shall adopt the relation for our notation, e.g., the 
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frequency at 92.908 cd" 1 has n — 32 and is called 

S32- 

There are several identical S sequences: the se- 
quences are also present as side lobes of the dom- 
inant modes. As an example, the dominant mode 
at fi6 = 66.297 cd _1 shows side lobes at 26.610 (ac- 
tually -26.610) and 159.206 cd _1 , corresponding 
to fi6-S32 and fi6+S32, respectively. There are 
a large number of additional agreements for all 
the dominant modes and S25 to S35. We empha- 
size that the large number of exact agreements 
between the observed and computed frequencies 
rules out accidental agreements and subsequent 
overinterpretation. 

The S sequences are particularly strong during 
Q2, in which fi also has a high amplitude. In sub- 
sequent quarters, the amplitudes of both fi and 
the S sequences decrease. This confirms the as- 
trophysical connection between the low and high 
frequencies. 

5.3. Triplet 3 and combinations with {4 

There exists a third connection between a dom- 
inant low frequency and the spacings of high 
frequencies. The triplet T3 between 47 and 
53 cd -1 has frequency separations of exactly 
(2.981 cd -1 ). Furthermore, various low-amplitude 
modes also show frequency differences correspond- 
ing to exactly {4. 

6. Spectroscopic observations 

Three spectra of KIC 8054146 were obtained on 
2011 June 29, July 16 and September 6 with the 
High Resolution Spectrograph (HRS) attached to 
the 9.2-m Hobby-Eberly Telescope (HET) at the 
McDonald observatory. This is a cross-dispersed 
echelle spectrograph yielding a resolving power 
(R) of 60000. The signal-to-noise ratio (S/N) per 
pixel at A ^5000A is 70, 82, and 75, respectively. 

The basic data reduction was carried out using 
the IRAjH software suite. Bias subtraction was 
performed by removing an averaged zero frame for 
each night, and spectral images of a flat- field lamp 
taken with the same setup as used for our target 



observations were used to correct uneven pixel re- 
sponse across the two CCDs. Wavelength calibra- 
tion (accurate to <~ 6 x 10~ 5 A) was derived from 
a ThAr emission lamp observed shortly before or 
after our stellar spectra were taken. Cosmic rays 
were identified as outliers to 2D fits to the spec- 
tral orders and removed. In rare cases where a 
cosmic ray fell across the echelle order, defeating 
our normal rejection routine, we have removed it 
with IRAF's COSMICRAYS utility. The finalized 
spectrum for each night was produced by coadding 
the three cosmic ray splits, resulting in a total ef- 
fective exposure time of 60 minutes. 

After the application of the heliocentric cor- 
rection, we checked that no radial velocity vari- 
ation was present among the three acquired spec- 
tra. This excludes the possibility that the star is a 
member of a close binary system. This allowed us 
to average the three acquired spectra to increase 
the S/N, reaching a S/N of 131 for the final spec- 
trum. 

The average spectrum, normalized by fitting a 
low order polynomial to carefully selected contin- 
uum points, covers the wavelength range 4100- 
7900 A, with a gap between 5927 A and 6024 A, 
because one echelle order is lost in the gap be- 
tween the two chips of the CCD mosaic detector. 
Given the very high projected rotational veloc- 
ity (wsini) of the star, we adopted the hydrogen 
lines as primary indicators for the effective tem- 
perature T e s determination. The normalization 
of hydrogen lines, observed with echelle spectro- 
graphs, gives the largest contribution to the total 
error budget in the T e g determination by hydrogen 
line fitting. To reduce this uncertainty we normal- 
ized the Ha and H/3 lines individually for each of 
the three spectra, obtaining in this way six quasi- 
independent measurements of T c ff. We normal- 
ized the hydrogen lines making use of the artificial 



flat- fie lding technique described in iBarklem et al 
(|2002h . 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



To compute model atmospheres, we employed 
the LLmodels stel lar model atmosphere code 
(jShulvak et al ] l2004h . For all the calculations, 
local thcrmodynamical equilibrium (LTE) and 
plane-parallel geom etry were assumed. We used 
the VALD database (Piskunov et al.lll995llKupka et al 



19991; iRvabchikova et al.l 119991 ) as a source of 
atomic line parameters for opacity calculations 
with the LLmodels code. Finally, convection was 



10 



i mplem ented according to the lCanuto fe Mazzitelli 
(Il991allbh model of convection. For more details 



sec 



Heiter et al 



42002D. 

We performed the T e g determination by fitting 
synthetic line prof iles, calculated with SYNTH3 
(jKochukhov 2007 ). to the observed Ha and H/3 
line profiles. For each spectrum we derived that 
the best fitting T c g is 7600 K for the Ha line and 
7700 K for the H/3 line, with a typical uncertainty 
of 150 K. 

To decrease even more the uncertainty on the 
effective temperature due to the normalization, 
on 2011 September 2 we obtained a spectrum of 
KIC 8054146 centered on the Ha line with the 
Cassegrain low-resolution spectrograph attached 
to the 1.8-m telescope of the Dominion Astrophys- 
ical Observatory (DAO), Canada. The adopted 
configuration of the spectrograph provided a re- 
solving power of about 15 000 and the obtained 
spectrum, covering the 6480 - 6780 A wavelength 
region, has a S/N per pixel of 60, calculated at 
about 6600 A. Since the spectrum covers the en- 
tire Ha line profile within a single order, the nor- 
malization line is a simple first order polynomial 
which can be safely determined by using the avail- 
able continuum regions on either side of the hydro- 
gen line. This decreases greatly the uncertainty 
on the measurement of T c g due to the continuum 
normalization. By fitting synthetic line profiles to 
this low-resolution spectrum, we derived a T e g- of 
7500 ± 200 K, where the uncertainty is mostly due 
to the quality of the spectrum. 

By averaging all these T e g measurements, we 
finally adopted T e g = 7600 K with a formal un- 
certainty of 80 K, which is not realistic, as de- 
scribed later. Figure [7] shows a comparison be- 
tween the observed HRS and DAO spectra with 
synthetic spectra, calculated with the adopted 
T cff of 7600 K and a T cS of 7400 K, for compar- 
ison. In the procedure of fitting synthetic hy- 
drogen line profiles to the observed spectra, we 
adopted a fixed value of log g = 3.9 (cgs) for the 
surface gravity and of 2.5kms _1 for the microtur- 
bulence velocity (w m ic)- All model atmospheres 
were pro duced assuming the s olar chemical com- 
position (jAsplund et all l2009h . With this set 



of fundamental parameters, by fitting the most 
prominent blends we produced synthetic spectra 
of the available spectral region and measured a 
radial velocity (u r ) of SliSkms" 1 and a vsini 




wavelength (A) 

Fig. 7. — Top panel: comparison between the Ha 
line profile (black solid line), observed with HRS, 
and synthetic profiles calculated with our best fit- 
ting effective temperature of 7600 K (red dashed 
line) and with an effective temperature of 7400 K 
(blue dotted line). For both synthetic spectra we 
adopted log g = 3.9, u m j C = 2.5kms _1 , and solar 
metallicity. Bottom panel: as in the top panel, 
but with the spectrum observed at the DAO ob- 
servatory. The sharp lines visible in the observed 
spectra are telluric spectral lines; the stellar pho- 
tospheric spectral lines are much shallower due to 
the high rotation. 



of 300 ±20 km s" 



Fol lowing the procedure described in lFossati et al.l 
(|2007l ). we attempted a measurement of the iron 
abundance from a few Fe I and Fe II lines, but the 
limited S/N of the spectrum, the extreme blending 
and the uncertainty on the continuum normaliza- 
tion prevented us from obtaining any useful result. 
As a consequence, we were not able to determine 
any reliable log g value by means of the ionisation 
equilibrium. At the temperature of KIC 8054146, 
the two analysed hydrogen lines display a little 
reaction to gravity variations, which allowed us to 
determine at least an upper and lower log g value 
of 4.2 and 3.4, respectively. The best fitting \ogg 
value is 3.9, indicating that the star is likely to be 
still on the main-sequence. A more precise \ogg 
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value could be obtained either with narrow band 
photometrjH, or with a distance value, or with cal- 
ibrated photometry spread over the region of the 
Balmcr jump, but none of this is currently avail- 
able. A more realistic value of ctt oH can then be 
obtained by taking into account the limits given 
for log g, leading to <TT cfl = 200 K. 

The large uncertainties on the abundance val- 
ues derived from the most prominent blends made 
it also impossible to determine a reliable v m i C 
value, which we therefore kept fixed at 2.5 kms" 1 



typic al of late A-type stars (see e.g. iFossati et al 



2008). Figure |8] shows a comparison between 



the LLmodels theoretical fluxes, calculated with 
T e g = 7600 K and log g = 3.9, and Johnson BV and 
2MASS JHK photometry. This comparison con- 
firms that the derived T e s fits well the available 
broad-band photometry and the absence of in- 
frared excess excludes the presence of warm dust 
around the star. 

In conclusion, the spectroscopic analysis al- 
lowed us to determine that KIC 8054146 is a very 
fast rotating late A-type star, which is likely to 
have a mctallicity close to that of the Sun. 

7. Discussion 

We have shown from the Kepler data that the 
low frequencies (up to 7cd _1 ) and the high fre- 
quencies (up to 200 cd" 1 ) are related in several 
families through preferred frequency spacings, fre- 
quency patterns, and similar amplitude variations 
in the individual families. Consequently, an expla- 
nation in terms of independent gravity and pres- 
sure modes appears too simple. The analyses of 
stars of spectral type A and F with the Kepler 
spacecraft have shown that the light o utput is 



modu lated by the rotational frequency ([Balona 



20111 ). Even a very 'simple' A star with radial 



pulsation modes and a very low rotation rate 
of ~ 20 kms" 1 already shows complex rotational 
variations including s mall amplitude mod ulations 



of the radial modes ( Breger et al.l 12011 ). Con 



scqucntly, a very rapidly rotating star such as 
KIC 8054146, presumably with differential rota- 



Uyttcrhoevcn ct al. (2011) published values of T e ff and 
log g obtained from Stromgrcn photometry, but a number 
of available narrowband measurements of this star have 
been found to be those of a nearby brighter star (TYC 
3149-1743-1). 




alength [A] 



Fig. 8. — Comparison between LLmodels theo- 
retical fluxes calculated with the fundamental pa- 
rameters derived for KIC 8054146 with Johnson 
BV and 2MASS JHK photometry. 



tion, might be very interesting. What is the rota- 
tion frequency of this star? 

Based on the mean values of T e g (7600 K) and 
log g (3.9) measured in the previous section we ob- 
tain a radius of about 1.9 Rq for the star (assum- 
ing a mass of 1.8 M Q ). Together with the mean 
measured value of the projected rotational veloc- 
ity, v sin i = 300 kms" 1 and assuming equator-on 
view (z = 90°) we derive a rotational frequency 
of 2.87 cd" 1 (corresponding to 0.75 of the Kep- 
lerian breakup rate). Taking into account the 
given uncertainty of 20 kms" 1 in vsini, we ob- 
tain a rotational frequency of 2.73 cd -1 (i.e., 0.69 
of Kcplcrian breakup rate) for the lower limit of 
v sin i — 280 kms -1 . The rotational frequency is, 
therefore, very similar to that of the dominant 
low frequencies between 2.8 and 3.0 cd" 1 , but the 
uncertainties exclude the possibility of identifying 
any one of the peaks as the rotational frequency. 

Due to the high value of the surface rotation 
velocity near breakup velocity, we can conclude 
that it is more likely that the star is viewed nearly 
equator-on. This fact helps us to reduce the un- 
certainties in the fundamental parameters, since, 
due to gravity darkening, parameters such as the 
effective temperature are a function of aspect in 
rapidly rotating stars. Consequently, the mea- 
sured temperature and gravity values are esscn- 
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tially the equatorial values. At the poles, we ex- 
pect gravity brightening, since both the T e g and 
logg values are higher there. 

7.1. Rotationally split 1=1 triplets? 

We do not have a physical explanation for this 
perfect link between the low frequencies and the 
higher-frequency triplets, but can offer a specula- 
tion: the triplets are the counterparts of the ro- 
tationally split I = 1 triplets commonly observed 
in slowly rotating stars. KIC 8054146, however, 
is an extremely fast rotator for which the amount 
of asymmetry in the frequency separations is not 
known. An effect that could make the frequency 
separation between rotation ally split modes sym - 
metric is the 1:1:1 resonance ( Buchler et al ] |l997ft . 
Such explanations in terms of resonances must, 
however, be regarded as speculative at this time. 

It is important to point out that the small, but 
statistically significant, departure from equidis- 
tance between the components of the two triplets, 
Ti and T 2 , is almost impossible to discover in 
stars not measured as extensively as KIC 8054146. 
The separations arc 2.9300 ± 0.0001 cd" 1 and 
2.9339±0.0002cd _1 , respectively. Even more re- 
markable is the fact that both these separations 
are seen as distinct low-frequency peaks and as a 
sum of these two values. Only the excellent fre- 
quency resolution and high accuracy of the two- 
year data made it possible to detect this. 

Could the low frequencies simply be combina- 
tion frequencies of the t—\ triplets? We have 
previously shown from long-term ground-based 
data for the 5 Set star 44 Tau, that the am- 
plitude changes of the combination modes mir- 
ror the amplitude ch anges of their parent modes 
( Breger fc Lend 2008 ) in the sense that the prod- 
uct of the parent amplitudes can be used to pre- 
dict the amplitude variations of the combination 
frequency. We have already presented arguments 
that the measured high v sin i value implies that 
we see the star almost equator-on. Consequently, 
photometrically the £ = 1, m = mode, which has 
a surface node line along the equator, almost can- 
cels out. The prograde and retrograde compo- 
nents, however, are fully visible. Therefore, the 
amplitude of the central component could be much 
higher than observed, raising the possibility of vis- 
ible combination frequencies involving the central 
component. 



Consequently, the low frequencies f 2 and f3 
can be combination frequencies, rather than grav- 
ity modes. Such an interpretation is fully sup- 
ported by the related quarterly amplitude variabil- 
ity shown in Figure UJ Even the puzzling result of 
Q10 is explained, viz. the amplitude of (f 2 +f3) 
has increased, while the amplitudes of both f 2 and 
f3 have decreased. In the combination frequency 
model, f 2 and f3 are the differences between m = 
and m = l, -1, respectively. Since the amplitude 
of the m = mode has decreased to nearly zero, 
the combinations are weak too. However, the am- 
plitudes of m~l and -1 have increased, so that 
the amplitude of the difference (=f 2 +f 3 ) has also 
increased. So the opposite behavior of f 2 and f3 
vs. (f 2 +f3) is explained. 

Of course, this remarkable agreement does not 
prove that the low frequencies are combinations, 
only that they behave like combinations, usually 
interpreted as nonlinearities or resonant excitation 
in the star. 

7.2. Eclipses? 

Could the frequency regularities be caused by 
eclipses? In the Fourier domain, the nonsinusoidal 
nature of eclipse light curves lead to long sequences 
of frequency multiples of the orbital frequency. In 
such cases inspection of the actual light curves 
may be more instructive. We have examined the 
light curves and find no evidence for eclipses. This 
does not prove that eclipses do not exist, since they 
could be hidden in the pulsations. 

There are two further arguments against eclipses: 
we have already determined that the S sequence 
with the equidistant fi spacing does not consist 
of exact multiples of f 1 . The offset argues against 
a Fourier series of a nonsinusoidal light curve. 
Furthermore, the values of the observed low fre- 
quencies are compatible with the size of the stellar 
radius, i.e., a stellar component would have to be 
unrealistically close to the main star. 

7.3. Spots and surface inhomogeneities 

We have already seen that the low frequencies 
fi to f4 are within the range of the predicted rota- 
tional frequency of the star. In cooler stars, such 
as the Sun, spots together with differential rota- 
tion cause light variations with similar frequencies 
and varying amplitudes. This is also observed for 
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KIC 8054146. We also note that the low frequen- 
cies are constant, i.e., if there were any spots, they 
would not be moving in latitude. Spots can also 
lead to asymmetric light curves, i.e., the observed 
frequencies between 5.5 and 6.0 cd -1 , but fail in 
the details such as the peaks not being simple har- 
monics of the lower frequencies. Even worse, spots 
cannot explain the observed close relationship be- 
tween the low and high frequencies. 

7.4. High-degree nonradial gravity modes 

Since KIC 8054146 is a rapid rotator, we have 
to take into account that the frequencies in the 
corotating reference system of the star may be 
quite different to the observed frequencies. The 
transformation from the corotating system to the 
observer's system causes the frequency separa- 
tion between rotationally split components to 
shift by an additional m/ ro t factor. In case of 
KIC 8054146, with / rot close to 2.8 cd" 1 this shifts 
prograde modes of higher m to rather high fre- 
quencies. 

If high-degree prograde gravity modes are ex- 
cited in the star, the equidistant high-frequency 
S sequence may in reality be a low-frequency se- 
quence shifted to the observer's frame of reference. 
In a subsequent investigation, we will model and 
test the work hypothesis of prograde high-degrec 
nonradial gravity modes (i.e., Kelvin modes) caus- 
ing the S sequences. 

7.5. Asymptotic pulsation 

It is not expected that the equidistant S se- 
quence with frequencies between 76 and 99 cd -1 
is due to acoustic asymptotic pulsation, because 
the He II zone cannot drive modes at such high 
frequencies. Moreover, no magnetic field has been 
detected, making roAp-type driving in the He I/H 
partial ionization zone unlikely as well. Further- 
more, while the distribution of mode frequencies 
within the asymptotic regime is predicted to be 
quite regular, the models d o not indicate exact 
equidistance. Reese et al. ( 20081 ) examined the 



acoustic mode frequencies in rapidly rotating stars 
based on polytropic models. As can be seen in Fig. 
1 in their paper, even at high radial orders there 
is still a small variation of the frequency spac- 
ings. This deviation from exact equidistance will 
increase if a more realistic stellar model (includ- 



ing the sharp chemical discontinuity at the border 
of convection zones) is used. Although the de- 
parture from equidistance is small enough not to 
destroy the regular patterns, it should be clearly 
significant with the accuracy of the given data set. 
Consequently, the exactly equidistant frequency 
spacings of the S sequence is unlikely to be due 
to asymptotic acoustic pulsation. 

7.6. Conclusion 

KIC 8054146 was shown to be a rapidly rotat- 
ing star near the cool edge of the instability strip 
(T c ff = 7600 ± 200 K, log g = 3.9 ± 0.3). The pho- 
tometric Kepler data covering two years revealed 
a large range of frequencies covering 200 cd -1 . Of 
the 349 statistically significant detected frequen- 
cies, a large group of 117 frequencies have val- 
ues less than 7cd _1 . Most of these are concen- 
trated in four specific frequency bands. They do 
not show the equidistant period spacings predicted 
for gravity modes of successive radial order, n, 
which was reported in at least one other hybrid 
pulsator. The four dominant low frequencies in 
the 2.8 to 3.0 cd -1 (32 to 35 fiRz) show strong am- 
plitude variability with a time scale of months and 
years. 

The four dominant low frequencies also deter- 
mine the spacing of the higher frequencies in and 
beyond the S Set pressure-mode frequency domain. 
In fact, most of the higher frequencies belong to 
one of three families with spacings linked to a 
specific dominant low frequency. In the Fourier 
spectrum, these family regularities show up as 
triplets, high-frequency sequences with absolutely 
equidistant frequency spacings, side lobes (ampli- 
tude modulations) and other regularities in fre- 
quency spacings. Within each family the ampli- 
tude variations between the low and high frequen- 
cies are related. We conclude that the low frequen- 
cies (gravity modes, rotation) and observed high 
frequencies (mostly pressure modes) are physically 
connected. This unusual behavior may be related 
to the very rapid rotation of the star. 
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